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1 Systems and States 

This lecture on polymer dynamics primarily addresses the molecular motion in melts and solu- 
tions. Other states like rubber or glass will be refered to briefly. Therefore, let us start by sorting 
out the following states of macromolecular matter: 

A polymer melt is basically a liquid. However, in contrast to a simple liquid, its constituent 
molecules are flexible, and therefore their center-of-mass motion is usually accompanied by 
conformational changes. The dynamics is further complicated by the entanglement of neigh- 
bouring chain molecules. As a result, polymer melts are highly viscous and non-Newtonian 
liquids. The prime example for a naturally occuring polymer melt is latex, a milky fluid rich in 
cis-polyisoprene that is found in many plants, particularly abundant in the rubber tree. 

A rubber (or elastomer) is obtained from a polymer melt by a chemical process called 
vulcanization that creates permanent inter-chain links. These cross-links prevent long-ranged 
diffusion and viscous flow. On short time and length scales, however, the polymer segments 
can move as freely as in a melt, which explains the characteristic elasticity of the rubber state. 
Everyday examples: erasers, car tyres. An unvulcanised polymer melt is said to be in a rubber 
state if the entanglement is so strong that the chains do not flow on the time scale of observation. 

Plastic materials are either thermosetting or thermoplastic. 

A thermosetting plastic is obtained by irreversibly creating strong chemical inter-chain 
links. Examples: bakelite, duroplast, epoxy resin. 

Thermoplastic materials are obtained in a reversible manner by cooling a polymer melt. 
They can be either amorphous or semi-crystalline. 

An amorphous solid is also called a glass, provided that it can be reversibly transformed 
into a liquid. This transformation is called the glass transition; it occurs at a loosely defined 
temperature Tg. Examples: poly (methyl methacrylate) (PMMA, plexiglas), polystyrene, poly- 
carbonate. 

In semicrystalline polymers, chains run through ordered and disordered domains. The or- 
dered, crystalline domains exist up to a melting point T^. They provide strong links between 
different chains. In the temperature range between Tg and Tm, the disordered domains provide 
flexibility. Example: polyethylene (Tg ~ -100 ... - 70°C, T^ ~ 110 . . . 140°C), used for foils 
and bags. Most other thermoplastics are typically employed below Tg. Example: polyethylene 
terephthalate (PET, Tg ~ 75°C, T„i ~ 260°C), used for bottles and textile fibers. 

In a polymer solution the dilute macromolecules are always flexible. Therefore they be- 
have as in a pure polymer melt, except that interactions between segments are modified by the 
presence of a solvent. 

2 Brownian Motion 

Almost all of polymer dynamics can be described by means of classical mechanics. Quantum 
mechanics seems relevant only for the explanation of certain low temperature anomalies that 

are beyond the scope of this lecture. Since we aim at understanding macroscopically relevant 
average properties of a large number of molecules, the appropriate methods are provided by 
statistical mechanics. 

Some fundamental methods of classical statistical mechanics have been originally devel- 
opped for the description of Brownian motion [1]. In this introductory section, we will review 
two of them: the Langevin equation for the evaluation of particle trajectories, and the Smolu- 
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chowski equation for a higher- level modelling in terms of probability densities. Both methods 
are ultimately equivalent, but depending on the application there can be a huge computational 
advantage in using the one or the other. ^ 



2.1 The Langevin equation 

We consider a colloidal (mesoscopic) particle P suspended in a liquid L. Molecules of L 
frequently collide with P, thereby exerting a random force F{t). On average the impacts from 
opposite directions cancel, so that 

{F{t)) = 0. (1) 

For this reason it was long believed that L — P collisions cannot be responsible for the random 
motion of P. By noting the analogy with the fluctuations of good or bad luck in gambling, 
Smoluchowski (1906) showed that this argument is fallacious: After n collisions, the expecta- 
tion value of total momentum transfer is indeed 0, but for each single history of n collisions one 
can nevertheless expect with high probability that the momentum transfer has summed up to a 
value different from 0. 

This random walk argument was put on a firm base by Langevin (1908).^ Let us write down 
the Newtonian equation of motion for the center-of-mass coordinate r{t) of P: 

md^r = -Cdtr + F, (2) 

or for the velocity v = dtv: 

mdtv = -Cv + F. (3) 

If P is a spherical particle of radius a, and L has a shear viscosity then the friction coefficient 
is (Stokes 1851) [4, § 20] 

C = 67r77sa. (4) 

Eq. (3) can be easily integrated: 

v{t)^m-^ f dt'F(Oe-^*-*')/^-. (5) 

J —oo 

P may be said to have a memory of past collisions that decays with a relaxation time 

Tm-^-^- (6) 

To compute averages like (v^) we must now specify the expectation value of random force cor- 
relation. We assume that different Cartesian components of F are uncorrelated with each other, 
{FaFp) = for a 7^ /3, and that random forces acting at different times t, t' are uncorrelated. 
The second moment of F then takes the form 

{Fa{t)Ff,{t'))=ASa0S{t-t') (7) 

where the delta function is a convenient approximation for a memory function that extends over 
no more than the duration of one L — P collision. 



'The following is standard material, covered by many statistical mechanics textbooks and by a number of 
monographs. For Brownian motion, see e. g. [2], for the Langevin equation, [3]. 

^K. Razi Naqvi [arXiv:physics/052141vl] contests this received opinion, arguing that Langevin's "analysis is 
at best incomplete, and at worst a mere tautology." 
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By virtue of (1), the average velocity is (v) = 0. However, as argued above, for each single 
history of F one can expect that the integral in (5) results in a nonzero velocity. After a short 
calculation we find that the second moment has the time-independent value 

- (8) 
On the other hand, this moment is well known from the equipartition theorem 

^{v') - Ik^T, (9) 

so that we obtain the coefficient A = 2kBT( that was left unspecified in (7). 

In the next step, we can compute the mean squared displacement of P within a time span t, 
which we will abbreviate as 

{r\t)) :^ {[r{t) - r{0)f) . (10) 
We determine r{t) by explicit integration of (5), and after some calculation [3] we find 

{r^{t)) = 6Dt - QDTra (1 - e"*/^") (11) 

where we have introduced 

D := (12) 

From the long-time limit {r'^{t)) ~ QDi, we identify D as the dijfusion coefficient (Sutherland 
1905, Einstein 1905). In the opposite limit t <^ r^, ballistic motion is found: (r^(t)) = 
{{v{0)tf). 



2.2 The Smoluchowski equation 

Alternatively, Brownian motion can be analysed in terms of the space-time distribution p(r, t) 
of suspended particles. The continuity equation is 

dtp + Vj^O. (13) 

As usual, the current j has a diffusive component (Fick 1855) 

jdiff ^ _£)Vp. (14) 

We now assume that a potential U (r) is acting upon the particles. In the stationary state, the 
force — VC/ is just cancelled by the friction —(v. Accordingly, the drift component of j is 

jdrift ^p^^ -C'pVU. (15) 

Collecting everything, we obtain the Smoluchowski equation (1915) 

dtp = DV^p + C'V{pVU). (16) 

For a velocity distribution or for the more general case of a phase-space distribution it is known 
as the Fokker-Planck equation (1914/1917). 
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Fig. 1: Temperature dependence of the specific volume of a glass forming system. Depending 
on the cooling rate, the cross-over from liquid-like to glass-like slopes is observed at different 
temperatures, and it results in different glass states. If the system is cooled rapidly to T2 and then 
kept at that temperature, it relaxes towards states that are normally obtained by slow cooling. 
Similar behavior is found for the enthalpy. 

In the stationary state dtp = 0, we expect to find a Boltzmann distribution p oc (— [///cbT). 
It is easily verified that this holds provided D and are related by the Sutherland-Einstein 
relation (12). Combined with (4), the Stokes -Einstein relation 



is obtained. Surprisingly, this relation holds not only for mesoscopic suspended particles, but 
even for the motion of one liquid molecule among others (the hydrodynamic radius a is then 
treated as an adjustable parameter). Only in the highly viscous supercooled state, the propor- 
tionality D (xT/r]s is found to break down. 

3 Segmental Relaxation and the Glass Transition 
3.1 The glass transition 

It is generally believed that the ground state of condensed matter is always an ordered one. 
However, in many supercooled liquids this ordered state is kinetically inaccessible. Crystalliza- 
tion through homogeneous nucleation requires the spontaneous formation of nuclei of a certain 
critical size. If the crystalline unit cell is complicated and the energy gain small, then this crit- 
ical size is rather larger. If at the same time a high viscosity entails a low molecular mobility, 
then the formation of a stable nucleus may be so improbable that crystallization is practically 
excluded. On further supercooling the liquid becomes an amorphous solid, a glass. 

Glass formation is observed in quite different classes of materials: in covalent networks 
like quartz glass (Si02) or industrial glass (Si02 with additives), in molecular liquids, in ionic 
mixtures, in aqueous solutions, and others. In polymers, glass formation is the rule, not the 
exception. As said in Sect. 1, some plastic materials are completely amorphous; in others 
crystalline domains are intercalated between amorphous regions. 

The glass transition is not a usual thermodynamic phase transition of first or second order, 
as can be seen from the temperature dependence of enthalpy or density (Fig. 1): the cross-over 



D 



(17) 
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Fig. 2: Dielectric spectrum in the viscous liquid or rubber state. As long as molecular dipoles 
are able to follow an external electric field modulation, the permittivity e(w) = e'(w) + ie"{uj) 
has the elevated low-frequency value of a polar liquid, here eo = 60. On the other hand, 
when probed with a high-frequency field modulation, dipoles appear frozen, resulting in the low 
permittivity of a solid, here eoo = 3. As required by the Kramers -Kronig relation, the dispersion 
of e'ioj) is accompanied by a dissipation peak in e"{uj). The dispersion step and the loss peak 
are stretched when compared to Debye's Eq. (82) (dashed). 



from liquid-like to solid-like slopes is smooth, and it depends on the cooling rate. When the 
cooling is interrupted within the transition range, structural relaxation towards a denser state is 
observed. This shows that glass is not an equilibrium state: on the contrary, it can be described 
as a liquid that has fallen out of equilibrium.^ 



3.2 Structural relaxation 

Relaxation can also be probed by applying a periodic perturbation. This is done for instance 
in dielectric spectroscopy (Fig. 2). For reference, the most elementary theory of dielectric 
relaxation (Debye 1913) is derived in Appendix B. It leads to a complex permittivity 

with a = 7 = 1. Empirically the cross-over from liquid-like to solid-like behavior extends over 
a much wider frequency range than this formula predicts: it is stretched. Often, this stretching 
is well described if the exponents anticipated in Eq. (18) are allowed to take values below 1 
(Cole, Cole 1941, Cole, Davidson 1951, Havriliak, Negami 1967). Alternatively (Williams, 
Watts 1970), e(ci;) can be fitted by the Fourier transform of the stretched exponential function 

<l>K(t)=exp(-(t/r)^) (19) 

originally introduced to describe relaxation in the time domain (R. Kohlrausch 1854, F. Kohl- 
rausch 1863). 



^You may object that the supercooled liquid is already out of equilibrium — with respect to crystallization. 
However, the crystalline ground state is a null set in phase space: beyond nucleation it has no impact upon the 
dynamics of the liquid. The glass transition, in contrast, can be described as an ergodicity breaking in phase space. 

The thermodynamics of the glass transition is subject of ongoing debate. A recent attempt to clarify confused 
concepts is [5]. 
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Fig. 3: (a) Correlation function t) of deuterated polybutadiene, measured by neutron spin- 
echo at q = 1.5 (close to the maximum of S{q) ). Solid lines: fits with fg exp(— (t/ r)'^) with 
fixed (3 = 0.45. (b) Master curve. Data from [6]. 

Generally, the relaxation time r depends much stronger on temperature than the stretching 
exponents 0,7 or (3. This can be expressed as a scaling law, often called time-temperature 
superposition principle: Permittivities e{uj; T) measured at different temperatures T can be 
rescaled with times r(T) to fall onto a common master curve e{T{T)uj). 

Quite similar results are obtained for other response functions. For instance, the longitudinal 
mechanic modulus, a linear combination of the more fundamental shear and bulk moduli, can be 
probed by ultrasonic propagation and attenuation in a kHz. . . MHz range, or by light scattering 
in a "hypersonic" GHz range (Brillouin 1922, Mandelstam 1926). 

Particularly interesting is the pure shear modulus G{co), which can be probed by torsional 
spectroscopy. Since a liquid cannot sustain static shear, there is no constant term in the low- 
frequency expansion 

G{uj) = ir]u + 0{u^). (20) 

The coefficient r] is the macroscopic shear viscosity. If G{uj) obeys time-temperature superpo- 
sition, then r]{T) is proportional to a shear relaxation time T,ri{T). 

By virtue of the fluctuation-dissipation theorem (Appendix A), relaxations can also be stud- 
ied in equilibrium, via correlation functions that can be measured in inelastic scattering experi- 
ments. Fig. 3 shows the normalized correlation function 

■.= S{q,t)/S{q,0) (21) 

of a glass-forming polymer, measured in the time domain by neutron spin echo. Scaling is 
demonstrated by construction of a master curve, stretching by fits with a Kohlrausch function. 

3.3 Relaxation map and secondary relaxation 

At this point it is interesting to compare the outcome of different experimental methods. It is 
found that the susceptibility master curve is not universal; different spectroscopies generally 
yield different stretching exponents. Relaxation times may vary by factors 2 or more. How- 
ever, with good accuracy all relaxation times have the same temperature dependence. This is 
demonstrated in Fig. 4 for a particularly well studied molecular liquid. Upon cooling from 290 
to 190 K, the relaxation times increase in parallel by almost 12 decades. When a viscosity of 
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Fig. 4: Relaxation times in glycerol, determined by various spectroscopies. This plot can be 
read as a dynamic phase diagram: whether the material reacts like a liquid or a solid depends 
not only on temperature, but also on the time scale of observation. Adapted from [7]. 



10^^ Poise (10^^ Pa s) or a relaxation time of about 10'^ s is reached, relaxation falls out of equi- 
librium on the time scale of a human observer. As anticipated above, this is the glass transition. 
All relaxation modes that follow the temperature dependence of viscosity are conventionally 
called (X relaxation. 

The temperature dependence of relaxation times is often discussed with reference to the 
simplest physical model that comes to mind: thermally activated jumps over an energy bar- 
rier Ea, 

T = To exp — — (22) 

(van t'Hoff 1884, Arrhenius 1889). Therefore, data are often plotted as logr versus l/T or 
Tg/T so that a straight line is obtained if (22) holds. The oc relaxation in glass forming liquids, 
however, almost never follows (22); its trace in the Arrhenius plot is more or less concave. Many 
fitting formulae have been proposed; for most applications it is sufficient to extend (22) by just 
one more parameter, as in the Vogel-Fulcher-Tammann equation (for polymers also named after 
Williams, Landel, Ferry) 

r = roexp- — — . (23) 
I — lo 

The singularity Tq lies below Tg, and it is unclear whether it has any physical meaning. 

In many materials, there is more than just one relaxation process. If the additional, sec- 
ondary processes are faster than a relaxation, they are conventionally labelled (3, y, In rarer 

cases, a slower process is found, designated as oc' and tentatively explained by weak intermolec- 
ular associations. While y and higher order relaxations are always attributed to innermolecular 
or side-chain motion, the situation is less clear for |3 relaxation. There are theoretical arguments 
(Goldstein 1969) and experimental findings (Johari 1970) to support the belief that (3 relaxation 
is a universal property of glass-forming systems. In any case, all secondary relaxations are 
well described by the Arrhenius law (22). This implies that at some temperature they merge 
with the a relaxation. This is well confirmed, principally by dielectric spectroscopy, which has 
the advantage of a particularly broad bandwidth. Fig. 5 shows an example for this merger or 
decoupling of cc and (3 relaxation; such an Arrhenius plot is also called relaxation map. 
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Fig. 5: cx and |3 relaxation times of cross-linked polyurethane, obtained by differential 
mechanical-thermal analysis (DMTA), photon correlation spectroscopy (PCS), dielectric spec- 
troscopy and neutron scattering. Data from [8]. 



3.4 The mode-coupling crossover 

At present, there exists no satisfactory microscopic theory of relaxation near the glass transition. 
At moderate to low viscosities the situation is slightly better, since basic phenomena can be 
explained to some extent by a mode-coupling theory (Gotze et al. 1984-)."^ This theory attacks 
the microscopic dynamics at the level of the normalized density pair correlation function (21). 
An equation of motion is written in the form 

= + z/,$,(t) + nl^,{t) + Ql [ dr^git - r)m,{$(r)}, (24) 

which guarantees that subsequent approximations do not violate conservation laws. In a system- 
atic, though uncontrolled expansion the memory kernel nig is then projected back onto products 
of pair correlations, 

m,{$(t)} ~ Yl Vkp,{T)Mt)%{t). (25) 

k+p=q 

The coupling coeffients Vkpq depend only on the static structure factor S{q), which in turn 
depends weakly on the temperature. This temperature dependence, however, is sufficient to 
trigger a transition from ergodic, liquid-like solutions to non-ergodic, glassy ones: 

$,(t^oo)^0 forT>rc, 
->cx)) -> /g > for T < Tc. 

On cooling towards Tc, a critical slowing down is predicted that is characterized by two diverg- 
ing timescales, 

° ' (21) 

with the reduced temperature a := T/Tc — 1. The microscopic timescale to is of the order 
The exponents fulfill < a < 6 < 1; they depend on just one lineshape parameter called A. 



"^The standard reference for mode-coupling theory is the comprehensive book [9]. More accessible introduc- 
tions are provided by [10, 11]. 
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Fig. 6: Incoherent intermediate scattering function Sseii{(l,t) of the molecular glass former 
ortho-terphenyl, obtained by combining Fourier transformed and resolution corrected neutron 
scattering spectra from three different spectrometers. Fits with the mode-coupling scaling law 
of fast (3 relaxation. Adapted from from [12]. 

The pair correlation function passes through the following scaling regimes: 

$,(t) ~ /, + /i,(TV2(t/t^)-- for to < t < t^, 

^ /, - hy/^B{t/Uf for U<^t< (28) 

The regime delimited by the first two equations has been been given the unfortunate name "fast 
|3 relaxation" (the Johari-Goldstein process is then called "slow |3 relaxation" though it is faster 
than a relaxation). 

The t^" power law is a genuine theoretical prediction; it has been searched for in many 
scattering experiments, and it has actually shown up in a number of liquids. With the power 
law, the theory explains experimental facts known since long (v. Schweidler 1907). This power 
law is also compatible with the short-time limit of Kohlrausch's stretched exponential; it leads 
over to the cc relaxation master curve implied by the third equation of (28). 

Mode-coupling predictions have been confirmed with impressive accuracy by light scat- 
tering studies of the density driven glass transition in colloidal suspensions (van Megen et al. 
199 1-). On the other hand, in conventional glass formers mode coupling is not the full story. By 
numerically solving simplified versions of (24), it is possible to fit relaxational spectra of nor- 
mal or slightly supercooled liquids. On further supercooling, in favorable cases (such as shown 
in Fig. 6) the power law asymptotes of slow |3 relaxation appear, and extrapolations yield a 
consistent estimate of Tc. Typically, this Tc is located 15% to 20% above Tg. This implies that 
the ergodicity breaking of Eq. (26) does not explain the glass transition; nor does the power law 
(27) fit the divergence of viscosities or relaxation times near Tg. 

This leads to the view that Tc marks a crossover between two dynamic regimes: a mode- 
coupling regime at elevated temperature and low viscosity, and a "hopping" regime in highly 
viscous, deeply supercooled liquids. Experimental support for the significance of this crossover 
comes from a possible change in the functional form of 'qiT) or r(T), from the aforementioned 
breakdown of the Stokes-Einstein relation (17), and from the merger of a and slow |3 relaxation. 
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Fig. 7: Spring-bead model of a flexible polymer. Since one bead represents several monomeric 
units, there is no preferential bond angle at the level of this model. Note also that the forces 
represented by the springs are entropic. 



4 Dynamics of a Free Chain 
4.1 The spring-bead model 

While the short-time dynamics of a polymer melt is quite similar to that of any other viscous 
liquid, on longer time scales the chain structure of the polymer makes a decisive difference, 
imposing strong constraints on the segmental motion. To study the motion of an entire chain, 
we neglect all details of chemical structure. We approximate the polymer as a sequence of 
beads at positions Vnin = 1, . . . , A^), connected by A^ — 1 springs. Each bead represents several 
monomeric units so that there is no preferred bond angle (Kuhn 1934[?], Fig. 7). 

The time-averaged equilibrium configuration of the polymer is assumed to be given by a 
Gaussian distribution of bead-to-bead vectors r„ — r„_i, 

N 



P{r} oc exp [ -K ^ (r„ - r^-i f ) , (29) 

where the force constant 



n=2 



^ := (30) 

ensures an average squared spring length of 6^. 

The size of a polymer coil can be characterized by a mean squared radius. The most com- 
mon measures are the end-to-end distance Re, 

Rl:={{rN-r,f) = Nb' (31) 

and the gyration radius 

^1 ■■= - ^g)') ^ Nb^G, (32) 

n 

with the center of mass 

rG:=N-'J2^n (33) 

The expressions (31), (32) hold in polymer melts and in solutions. In other cases, both 
expressions are still valid approximations if the factor is replaced by N'^'^. In good solvents, 
the exponent is z/ = 3/5. 
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4.2 The Rouse model 

The Gaussian distribution (29) is based on the assumption that the free energy A = U — TS is 
dominated by the entropy S = kslnP so that the internal energy U can be neglected. Hence 
each bead experiences an entropic force 

FT' = --^A = -K i-Tn-l + 2r„ - Vn+l) , (34) 

with obvious modifications for n — 1, N. This force strives to minimize the distances between 
beads, thereby maximizing the coiling of the polymer. 

The coupling to the heat bath shall be modelled by a random force F^'^^. Its second moment 
is given by an obvious extension of (7), 

(F„^r'(*)0'(^')) = 2Ck^T5nmSa0Sit - t')- (35) 

Finally, moving beads experience a friction —C^dtTn. These three forces make up the Rouse 
model, which is a key reference in polymer physics (Rouse 1953).^ Accordingly, the Langevin 
equation is 

mdlrr, = -CdtVn + + F'^:-'. (36) 

In polymer solutions, the simple linear friction term is no longer adequate; it must be replaced 
by a hydrodynamic interaction. This interaction is handled reasonably well by a theory (Zimm 
1956) outlined in Appendix C. 

At this point, let us indicate some orders of magnitude. Assuming a spring length h = 
1 nm and equating it with the hydrodynamic radius in (4), and assuming further a microscopic 
viscosity % ~ 10 Pa-s, we find a friction coefficient ( of the order of 10~^ to 10~^ Ns/m, 
in agreement with empirical data for polyisobutylene, polymethyl acrylate and natural rubber 
[16]. Assuming a bead mass of 100 Da, the single-bead collision relaxation time (6) is = 
m/( ~ 10^^*^ s, which means that inertia is completely negligible on all relevant time scales. 
At T = 300 K, thermal motion is of the order (v^)^/^ ~ 300 m/s, and the force constant is about 
K ~ 10-2 N/m. 



4.3 Rouse modes 

We will solve (36) by transforming to normal coordinates. To begin, we note that there is no 
coupling whatsoever between the three Cartesian components. Therefore, we need to consider 
just one of them. Let us write x and / for an arbitrary component of r and F^^^^. Then the 
Langevin equation reads 

md^Xn = -(dtXn - K {-Xn-1 + 2Xn - Xn+l) + fn- (37) 

Introducing the vector notation x:—{xi,..., xn)'^, Eq. (37) takes the form 

md^x = —(dtx — kK_x + / (38) 

^Rouse theory is covered in many text books, most often using a continuum approximation probably due to de 
Gennes [13]. The short but well written chapter in [10] comes with a nice selection of experimental results. Among 
dedicated polymer physics textbooks, I found [14] indispensable though largely indigestible for its coverage of 
computational details, and [15] inspiring though sometimes suspicious for its cursory outline of theory. 
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NxN 

The entries +1 at both extremities of the diagonal reflect the necessary modification of Eq. (37) 
for n = 1,N. In the well known derivation of phonon dispersion, this complication at the 
boundary is usually ignored or superseded by an unphysical periodicity. It is largely unknown 
that the correct X can be diagonalized quite easily without any approximation. The eigenvalues 
are 



A. 



2 - 2 cos 



2 PTT 



4 sin - — , p 



and the normalized eigenvectors have components 

7V-V2 



0,...,N-1, 



forp = 0, 



(40) 



, . I/O pin — k)'^ 
(Ar/2)~^/^cos ^ for all other p, 



(41) 



where n = 1, . . . , A^. The proof requires no more than a straightforward verification of I^Vp = 
\pV. Collecting the normalized eigenvectors into an orthogonal matrix A :— (vq, . . . , ^ jv-i) 
we introduce normal coordinates 



X 



Ax, 



(42) 



and similar for /. It is easily seen that the average random force correlation (35) is still diagonal, 

if pit) flit')) = 2CkBTSpgSit - t'). (43) 
In consequence, for each normal mode one obtains a decoupled Langevin equation 

mdtXp = -(dtXp - nXpXp + fp. (44) 

At this point we must distinguish the eigenmode with the special eigenvalue Aq = from all 
the others. 

The eigenmode p — describes the motion of the center of mass vq — N'^^'^fo. Since 
Ao = 0, (44) is identical with the Langevin equation for Brownian motion (2). Accordingly, the 
long-time evolution of the mean squared displacement is (rQ(i)) ~ QD^tt, where the macro- 
molecular Rouse diffusion coefficient is given by a rescaled version of the Sutherland-Einstein 
relation (12): 

Turning to the Rouse modes with p > 0, we neglect the inertial term in the Langevin 
equation (44). Integration is then straightforward: 



Xpit) = C' dt'e-(*-*')/^-/,(0, 



(46) 
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Fig. 8: (a) Cis-polyisoprene consists of monomers that possess a dipole moment in direction 
of the polymeric bonds, (b) Peak frequencies u of dielectric loss spectra in cis-polyisoprene 
samples of different molecular weight M. The normal mode resonance depends strongly on 
M, whereas the a relaxation does not. (c) Normal mode relaxation times t = u^^ versus M, 
showing a crossover from the Rouse regime ~ to a stronger M dependence above the 
entanglement onset M^. Data from [17]. 

introducing the mode relaxation time Tp := C/{kXp). It can be approximated fox p <^ N as 

Tr, ^ ^. (47) 
The relaxation time of the fundamental mode p = 1 is known as the Rouse time: 

TR := t^^tN^. (48) 

In this expression, and h enter only via the extended chain length L = Nb, which does not 
change if we change our spring-bead model to consist of N/x segments of length bx. This 
justifies ex post that we have ignored all details of microscopic structure. A simple alternative 
expression for tr is 

2Rl 

TR := (49) 

with the gyration radius (32). 



4.4 Macroscopic consequences: dielectric and shear response 

In favorable cases. Rouse relaxation can be observed by dielectric spectroscopy. In the simplest 
case, the monomeric possess a dielectric moment ^ in direction of the polymeric bond. This 
requires the absence of a perpendicular mirror plane, which is the case e. g. for cis-polyisoprene 
(Fig. 8a). Then the overall dipole moment of the macromolecule is 

N 

= ~ ■f'n-O/U = (rN - ri)/i (50) 

n=2 

From (41) we infer that only eigenmodes with odd p contribute to (50). For t > tr pair 
correlations are dominated by the p = 1 mode: 

(^(O)Ai(t)) oc (fi(O)fi(t)) cx e-*/^. (51) 
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According to the fluctuation-dissipation theorem (73), this correlation is proportional to a linear 
response function. After Fourier transform, one finds that the dielectric permittivity e{uj) has a 
Debye resonance around a; ~ r£^. 

Fig. 8b,c shows results of dielectric spectroscopy in cis-polyisoprene melts with different 
extended chain lengths L = Nb. In experimental reports the chain length is of course expressed 
as molecular weight M. The a. relaxation peak, discussed above in Sect. 3.2, is perfectly inde- 
pendent of M, which confirms that it is due to innersegmental motion. In contrast, the normal 
mode relaxation time depends strongly on M. For low M, the Rouse prediction tr ~ is 
confirmed. However, at Mc ~ 10 kDalton, there is a rather sharp crossover to a steeper power 
law M^^ that is ascribed to entanglement effects. 

Mechanical spectroscopy has the advantage that it works also if monomers are too sym- 
metric for dielectric measurements. As already mentioned in Sect. 3.2 torsional spectroscopy 
probes a system's response to shear. Not seldom this response is nonlinear, due to non- 
Newtonian flow phenomena that are beyond the scope of the present lecture. As long as the 
response is linear, it can be described by the frequency dependent shear modulus G{u). Its high 
frequency limit Goo quantifies the shear needed to cause a given strain. In a liquid, the low fre- 
quency limit Go is zero: stationary shear is not able to build up a lasting stress. Sometimes this 
is seen as the defining property of the liquid state. Instead of stationary shear, a flow gradient is 
needed to maintain strain. The proportionality coefficient is the shear viscosity rj, which is the 
low frequency limit of G{uj) / (ico), as anticipated in (20). 

To calculate G{u) in the frame of Rouse theory, it is convenient to invoke a Green-Kubo re- 
lation (a variant of the fluctuation-dissipation theorem) according to which G{uj) is proportional 
to the Fourier transform of a stress autocorrelation function {(Jxz{t)<^xz{^)) The stress compo- 
nent Gxz can be expressed through the displacements of individual beads. For weak elongations, 
the autocorrelation can be factorized so that G{t) is proportional to the square of a normalized 
one-dimensional displacement autocorrelation function 

^ {5Xn{t)5Xn{Q)) ^ ^ {Xp{t)Xp{0)) 

Using (46) to compute the normal mode autocorrelation we find 

G(t) - J]e-'*/^-. (53) 
p 

With Tp oc and replacing the sum by an integral we obtain the approximation 



G{t) ~ Ve-^f'*/^ ~ / ~ ( 

P>i 

Fourier transform yields the power law 



(54) 



G{uj) - o;^/^. (55) 

For small t, quite many eigenmodes contribute to (53) so that neither the low-p expansion (47) 
nor the extension of the summation to oo are justified. Therefore G{uj) must cross over from 
(55) to a constant high-frequency limit Goo- 



^The following is no more than a speculative summary of obscure calculations in [14] and [15]. 



IFF Spring School 2011, B4. 16 



Joachim Wuttke 




10 



1-3 



10° 

(0 (s"^ at 1 60 K) 



Fig. 9: Real part G'{uj) of the shear modulus of polystyrene samples with different narrow 
molecular weight distributions. Experiments were done in a frequency range from 10^^'^ to 
j^go.s ^-1 ^-^ ^ temperature range from 120 to 260 K. Then, time-temperature superposition 
was used to construct the master curves shown here. Data from [18]. 

For large t, only few eigenmodes contribute to (53) so that the passage to a continuous p 
becomes invalid. In this case, it is more appropriate to compute the Fourier transform term by 
term, which yields a sum of Maxwell-Debye resonances 



In the low-frequency limit we obtain a constant term of doubtful significance plus a linear term 
proportional to the viscosity 



From (48), we have tr ~ iV^, but there is a prefactor N^^ in G{uj) omitted in our sloppy 
derivation so that finally the Rouse model predicts 77 ~ A^. 

In Fig. 9 experimental data are shown. For moderate chain lengths, G'{uj) is indeed found 
to cross over from co"^ (the lowest non-constant real term in the expansion of (56) to w^/^. The 
ultimate limit Goo has not been reached in this experiment. For longer polymer chains, a flat 
plateau appears between the liquid-like and the Rouse regime w^/^. Such a constant value of 
G{uj) implies instantaneous, memory-free response, which is characteristic of rubber elasticity; 
it is caused by entanglement. 

4.5 Microscopic verification: neutron spin echo 

For direct, microscopic measurement of chain conformation fluctuations one must access length 
and time scales of the order of nm and ns. The most powerful instrument in this domain is the 
neutron spin echo spectrometer. The recent book [19] provides a comprehensive review of spin 
echo studies on polymer dynamics. 

Usually, spin echo experiments require deuterated samples to avoid the otherwise dominant 
incoherent scattering from protons. The experiments then yield the coherent dynamic structure 
factor S{q,t). However, for a simple, intuitive data analysis the incoherent scattering func- 
tion S\{q, t) is preferable. It is also denoted Sscif (g, t) since it reveals the self correlation of a 




(56) 




(57) 
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Fig. 10: Single-chain coherent normalized dynamic structure factor of polyethylene melts at 
509 K. Lines are predictions of the Rouse model. They fit for short chains (left), but not for long 
chains (right); note the different time scales. Data points and lines from [19]. 

tagged particle. In Gaussian approximation 

SM.t) = (exp {iq(r{t) - r(0)))) - exp {-q^{r\t)) /Q) , (58) 

it yields the mean squared displacement (10). 

Measuring Si{q,t) by neutron spin echo is difficult because the random spin flips associ- 
ated with incoherent scattering destroy 2/3 of the incoming polarization. Nevertheless, thanks 
to progress in instrumentation, it is nowadays possible to obtain self correlation functions from 
undeuterated ("protonated") samples in decent quality. Alternatively, self correlations can be 
measured with the data quality of coherent scattering if short protonated sequences are interca- 
lated at random in deuterated chains. 

Within the Rouse model, and neglecting ballistic short-time terms, the mean squared dis- 
placement is given by 

AT-l C N-1 ^ 

= ^ E - (^p(t)^pm] = 6Dr t + 5^ [1 - e-*/^^] . (59) 

p=0 I p=l J 

With the same technique as in (54), one obtains the approximation 

(r2(t))^6DR|t + (7rrRt)^/'}. (60) 

At about t ~ tr, there is a cross-over from a t^/^ regime dominated by conformational fluctu- 
ations to the diffusion limit. Inserting the asymptotic (r^(t)) ~ t^/^ into (58) one obtains an 
expression that agrees with the Kohlrausch function (19) with a stretching exponent (3 = 1/2. 
This indicates that the high-p limit of the Rouse modes is more physical than might have been 
expected; it seems to capture even some aspects of segmental a relaxation. Be that as it may, 
the t^/^ prediction has been impressively confirmed in neutron scattering experiments. 

The coherent dynamic structure factor is more involved than (58). In general, it contains 
contributions from interchain as well as from intrachain correlations. Single-chain dynamics 
can be isolated by contrast variation, using a mixture of about 10% protonated and 90% deuter- 
ated polymer. Fig. 10 shows the single-chain dynamic structure factor of two polyethylene 
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Fig. 11: Entanglement, (a) Local detail, with relatively few chain crossings, (b) Coarse grained 
representation, showing heavy entanglement. Each monochrome, tube-like region is filled with 
coiled subunits of the polymer For times between the entanglement time tg and the disentangle- 
ment time td, polymer motion is confined to these tubes. 

melts with different chain lengths [19]. For short chains, the data are in perfect agreement with 
the Rouse model. For long chains, however, correlations decay much slower than predicted by 
the Rouse model. This is yet another indication of entanglement. 

5 Entanglement and Reptation 

Entanglement means that the conformational dynamics of a chain is hindered by presence of 
other chains. Entanglement is a topological constraint, due to the simple fact that chains cannot 
cross each other (Fig. 11). We have already encountered experimental results that provide clear 
evidence for the relevance of entanglement for polymer chains that exceed a certain size Nc'. 

• The dielectric relaxation time crosses over from the Rouse behavior r oc A^^ (48) to a 
steeper slope r oc A^'^'^ (Fig. 8c). 

• In the shear modulus G{u), there appears a plateau G{co) = const with rubber-like elas- 
ticity between the liquid limit G{uj) ~ ir]u and the Rouse regime G{uj) ~ o;^/^ (Fig. 9). 

• Neutron scattering shows that correlations within long chains decay much slower than 
predicted by the Rouse model (Fig. 10). 

The rather sharp crossover at Nc implies that entanglement becomes relevant if the coil ra- 
dius N^/'^h (up to a constant factor, depending on definition) exceeds a certain value 



Up to this length scale, the chains are heavily coiled with little mutual penetration (Fig. 11a). 
On larger scales, the coarse-grained polymer chains have the character of heavily entangled 
tubes (Fig. 1 lb). Each tube can be modelled as an ideal random chain, consisting of beads of 
size a. 

An entanglement time can be defined by (r^(te)) — a^. Using (60) in the limit t <^ tr, 
we find up to numeric prefactors 



with the critical extended chain length Lc = Neb. For times beyond te, the dynamics is qualita- 
tively different from the free chain Rouse regime. Since a chain is basically confined to a tube, 
it can only perform a one-dimensional, snake-like motion, called reptation (de Gennes 1971). 



a : = 



N'J^b. 



(61) 




(62) 
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Fig. 12: Time evolution of the mean squared displacement (10) on a double logarithmic scale, 
(a) For short chains, as predicted by Rouse theory, (b) For long chains, as predicted by de 
Gennes ' reptation theory. 



For a short outline of some scaling results, we concentrate on the mean squared displace- 
ment. We will see that there are altogether no less than five different regimes. They are summa- 
rized in Fig. 12. 

The one-dimensional dynamics within a tube shall be described by the Rouse model as 
before. Let s be a coordinate along the tube. The mean squared displacement in s is just one 
third of the Rouse result (60) in r. Since the tubes are ideal Gaussian random coils, an extended 
tube length = N^a^ corresponds to a squared real-space displacement of = NgO^ = as. 
In the reptation regime tc ^ ^ ^ we obtain, omitting prefactors, 

{r'{t)) ~ aDK\Tj,t f\ (63) 

This t^/^ law is a key prediction of reptation theory; it has been confirmed by neutron spin echo 
measurements [19]. 

For times beyond t^, the one-dimensional dynamics crosses over from innerchain Rouse 
fluctuations to center-of-mass diffusion. Accordingly, the real-space mean squared displace- 
ment takes the form 

{r\t)) ~ aDllY/'. (64) 

This holds until the chain escapes from its tube, which happens when ~ N^b'^ or ~ 
aD^t]/^. Using again = as, we obtain the disentanglement time 

td ~ (65) 

Finally, on time scales above t^, the chain, having diffused out of its original tube, is free 
to try new conformations in three dimensions. This is a center-of-mass random walk, described 
by 

{r\t)) ^ 6Ddt. (66) 



Matching (66) with (64) at t^, we obtain the disentangled diffusion coefficient 
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The scaling laws ~ and ~ N~'^ are important predictions. The disentanglement 
time td determines the relaxation time r observed in dielectric or mechanical spectroscopy. Em- 
pirically, the molecular mass dependence of r is even stronger than N^. Typical exponents are 
3.2 to 3.6; in Fig. 8c we had even 3.7. This discrepancy shows that one-dimensional diffusion 
in fixed tubes is not the full story. It is necessary to take into account fluctuations of the neigh- 
bouring chains (contour length fluctuations). On the other hand, the prediction ~ has 
been confirmed by quite direct, spatially resolved diffusion measurements [15]. 



Appendices 

A Linear response theory: relaxation, dissipation, fluctua- 
tion 

If a multi-particle system is exposed to a weak perturbation A, its response B is linear in A, 
as far as amplitudes are concerned. However, the response may be delayed in time, assuming 
the character of relaxation. Relaxation may be probed in time or in frequency, by spectroscopy 
(response to external perturbation) or by scattering methods (fluctuations in equilibrium). The 
relations between these probes are the subject of linear response theory, to be briefly summa- 
rized in this appendix. 

The linear response B{t) to a perturbation A{t) can be written as 

B{t)= f dt' R{t-t')A{t'). (68) 

J — oo 

Consider first the momentary perturbation A{t) = 6{t). The response is B{t) = R{t). There- 
fore, the memory kernel R is identified as the response function. 

Consider next a perturbation A{t) — e''*©(— i) that is slowly switched on and suddenly 
switched off (9 is the Heavyside step function, r] is sent to 0+ at the end of the calculation). For 
t > 0, one obtains B{t) — $(i) where $ is the negative primitive of the response function, 

R{t) = -dMt). (69) 

Since $ describes the time evolution after an external perturbation has been switched off, it is 
called the relaxation function. In the special case of exponential (Lorentzian) relaxation, R and 
$ are equal (up to a constant factor), which is a frequent source of confusion. 

Consider finally aperiodic perturbation that is switched on adiabatically, A{t) = exp{—iujt+ 
rjt), implying again the limit r] — )■ 0+. The response can be written B{t) = x(c<;)yl(t), introduc- 
ing a dynamic susceptibility 

X{uj) := / dte'^'^^''^^* R{t). (70) 
Jo 

This motivates the definition 

poo 

F{uj) := / dte*'^*$(t). (71) 
Jo 

of the one-sided Fourier transform F{uj) of the relaxation function ^{t). Because of (69), there 
is a simple relation between x and F: 



XH = ^0)+iujF{uj). 



(72) 
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In consequence, the imaginary part of the susceptibility, which typically describes the loss peak 
in a spectroscopic experiment, is given by the real part of the Fourier transform of the relaxation 
function, Im % = wRe F{uj). Conversely, dispersion is described by Re % = $(0) — wim F{u). 

Up to this point, the only physical input has been Eq. (68). To make a connection with 
correlation functions, more substantial input is needed. Using the full apparatus of statistical 
mechanics (Poisson brackets, Liouville equation, Boltzmann distribution, Yvon's theorem), it 
is found [20] that for classical systems 

{A{t)B{Q))^k^T^{t). (73) 

This is an expression of t\\Q fluctuation-dissipation theorem (Nyquist 1928, Callen, Welton 
1951): the left side describes fluctuations in equilibrium; the right side relaxation towards equi- 
librium, which is inevitably accompanied by dissipation (loss peak in Im x). 

Pair correlation functions are typically measured in scattering experiments. For instance, 
inelastic neutron scattering at wavenumber q measures the scattering law uj), which is the 
Fourier transform of the density correlation function, 

1 r°° 

S{q,uj) = —J dte^-'{p{q,trp{q,0)). (74) 

In contrast to (71) and (70), this is a normal, two-sided Fourier transform. In consequence, if 
we let {p{q, t)*p{q, 0)) = ^{t), then the scattering law S{q, u) is proportional to the real part 
Re F{u!) of the one-sided Fourier transform of ^{t). 



B Debye's theory of dielectric relaxation 

In modem terms, Debye's theory of dipolar relaxation is based on a Smoluchowski equation 
that describes the time evolution of the probability distribution /(^?, 0, t) of dipole orientations 

^,(j)(0 <^ <7i,0 <(f)< 2n): 

C9t/ = V(rV/ + C/V/) (75) 

where C is a friction coefficient, U is an external potential, and /3 = 1/ {k-oT). Inserting spherical 
coordinates, ignoring 0, and keeping r constant, we obtain 

C'dtf ^ {sm^y^d^sin^ir^d^f + fd^U) (76) 

with C,' = r'^C,. An electric field in z direction causes a potential 

U(t) ^-pE(t) cos '& (77) 

that is proportional to the dipole moment p. 
With the ansatz 

mt)^l + (3g{t)cos^, (78) 
and introducing the relaxation time r := /3CV2, Eq. (76) simplifies to 



rdtg{t) = 



-g{t)+pE{t) + 0{gPpE). 



(79) 
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Under realistic experimental conditions, we always have fiE <^ (3^^ so that the last term is 
negligible. The remaining linear differential equation shall be rewritten for a macroscopic ob- 
servable, the polarization 

P(t)^j^,oosm«,t) = ^- (80) 

We obtain 

{l + rdt)P{t)^^E{t). (81) 

In the simplest time-dependent experiment, the electric field is adiabatically switched on, 
then suddenly switched off at t = 0. The polarization then relaxes exponentially, P{t) oc 
exp(— t/r). In a frequency-dependent experiment, a periodic perturbation E{t) oc exp{iu!t) is 
applied. This yields the susceptibility 

Xdipolar(w) = = . (82) 

E{UJ) 1 — lOUT 

The relative electric permittivity is then 

e{uj) = 1 + Xdipolar(w) + Xothcr(w) (83) 

where the "other" contribution comes mainly from the electronic polarizability. The dashed 
lines in Fig. 2 show the dispersion step in the real part e'{uj) and the dissipation maximum in 
the imaginary part e"{u). 



C Zimm's theory of chain dynamics in solution 

Starting with equation (34), the entire Rouse model is based on the assumption that the chain 
conformation is driven by entropy only. This a good approximation for melts, but generally not 
for solutions, except in the condition. To account for the swelling of a polymer in solution, 
one needs at least to model the mutual sterical exclusion of different chain segments. The 
simplest approximation for this excluded volume interaction is the repulsive potential 

). (84) 

As described in Sect. 4.1, its effect upon the equilibrium structure is limited to the modification 
of the exponent u in the scaling laws (31), (32) for the coil radius. 

For the dynamics, another modification of the Rouse model is even more important: one 
has to include the hydrodynamic interaction between the polymer and the solvent. The motion 
of a polymer bead drags the surrounding solvent with it, thereby creating a flow pattern, which 
in turn exerts a force upon other beads. If inertia is neglected, the friction term assumed in the 
Rouse model implies v — (~^F.To account for hydrodynamic interactions, this equation must 
be replaced by 

Vn = ^ HnmFm- (85) 

m 

To estimate the coupling coefficients Hnm, one usually refers to a simple case for which the 
hydrodynamic interaction can be obtained from first principles: for a point particle, located at 
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ri{t) and dragged by a force Fi, one can solve the Navier-Stokes equations to obtain the flow 
field v{r) = H{r — ri)Fi with the Oseen tensor^ 

H{r) = (1 + f ® f) , (86) 

which is then used to approximate 

f°^^ = ^' (87) 

\ H{rn - Tm) else. 

Unfortunately, the r dependence of H makes the modified Langevin equation nonlinear. This 
obstacle is overcome in the Zimm theory (1956) by a preaveraging step that is basically a mean 
field approximation: Hnm is replaced by its average under the equilibrium distribution P{r}. 
In the © condition, one obtains forn 7^ m 

(Hnm) ^ 3 ^ , (88) 

(D7r'*|n — m\) ' rjo 

which is a rather long-ranged interaction. In other conditions, a modified distribution P{r} 
might be used, leading to a modified power law \n — m\^'^. 

The preaveraging linearization allows to rewrite the one-dimensional Rouse mode Langevin 
equation (44) with hydrodynamic interaction as 

dtx = K (y-mdlx -KM + lj (89) 

with the diagonal matrix A^^ — 5pqXp and with — I^A, which in the 6 condition is in 
good approximation 

'0 for p 7^ g, 
-J- — forp = g = 0, 

2>{<o7i^f'^r]h (90) 
iVV2 



H ~ < 
PI — » 



else. 



For the p = eigenmode, we find again Brownian motion, with the Zimm diffusion constant 

= oc N-^'\ (91) 

which differs from Dr ~ in the Rouse model. Allowing for excluded volume interaction in 
a good solvent, the aforementioned generalisation of (88) leads to Dz ~ N~^. Comparing with 
(32), we find that the diffusion constant is in both cases determined by the coil radius: Dz ~ 
R^^. This result is routinely used in photon correlation spectroscopy (somewhat misleadingly 
also called dynamic light scattering), where the diffusion coefficient of dilute macromolecules 
is measured in order to determine their gyration radius. 



am unable to trace back this result to C. W. Oseen whose 1910 papers in Ark. Mat. Astr. Fys. are frankly 
unreadable. Zimm (1956) takes the tensor from Kirkwood and Riseman (1948) who cite a report by Burgers (1938) 
that is not easily available. 
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For nonzero eigenmodes, the Rouse time is replaced in the 6 condition by 

-~ (37r)V2A;BT' ^ ^ 

and the mode relaxation times become 

- (93) 
Again, the N dependence can be generalized towards a dependence on the coil radius, tz ~ -R^. 
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